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CHAPTER  I 


INTRODUCTION 

The  transmission  of  infrared  radiation  through  the  atmosphere  is 
of  considerable  interest  to  many  disciplines  of  science  and  a  drive  to 
better  understand  the  various  attenuation  processes  of  molecular 
absorption  and  aerosol  scattering  constantly  goes  on.  These  attenuation 
processes  are  particularly  dependent  on  tropospheric  concentrations  of 
Hod  and  CO2  as  well  as  aerosol  concentrations.  Water  vapor  absorption 
dominates  much  of  the  infrared  spectrum  and  affects  the  performance  of 
many  E-0  systems.  This  absorption  is  directly  dependent  on  the 
concentration  of  water  vapor  molecules  along  the  atmospheric  path  of 
interest. 

Some  of  these  E-0  systems  include  laser  radar  or  LIOAR  systems 
which  can  measure  atmospheric  water  vapor  profiles,  temperature  profiles 
and  gaseous  pollutant  concentrations.  Other  infrared  systems  include 
optical  data  transmission  links,  laser  range  finders,  and  infrared  laser 
target  designators.  Some  passive  infrared  systems  include  thermal 
imagers  and  Forward-Looking  Infrared  (FUR). 
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These  systems  are  largely  confined  to  the  optical  windows  or 
regions  where  the  attenuation  is  small.  Figure  1.1  shows  that  the 
infrared  windows  of  the  troposphere  occur  between  IT  to  8  um,  5  to 
4.5  um  and  4.1  to  3.2  um,  and  are  bounded  by  strong  HjO  and  CO? 
absorption  regions.  The  weak  absorption  of  the  window  regions  is  caused 
by  local  line  and  continuum  absorption. 

Continuum  absorption  is  a  depression  of  the  spectral  background 
throughout  the  infrared  spectrum  which  appears  as  a  slowly  varying 
function  of  frequency.  While  many  continuum  studies  have  been  done  in 
the  10  um  region  [1,2, 3, 4, 5],  not  nearly  as  much  has  been  done  in  the  3 
to  5  um  band  [6,7].  This  is  mostly  because  the  continuum  absorption  in 
this  region  is  an  order  of  magnitude  smaller  than  the  continuum 
absorption  in  the  8  to  13  um  region.  Although  there  is  greater 
continuum  interest  in  the  8  to  13  um  band,  more  data  is  needed  from  the 
3  to  5  um  band  to  gain  a  clearer  picture  of  the  entire  continuum  effect. 

In  1938,  based  on  his  observations  of  the  8  to  12  um  band, 

W.M.  Hlsasser  [8]  suggested  the  existence  of  a  water  vapor  continuum 
and  speculated  that  it  was  caused  by  the  far  wings  of  very  strong  water 
vapor  lines  outside  the  window  region.  However,  a  line  by  line 
computation  of  the  continuum  absorption  by  Burch  [2],  using  a  Lorentzian 
line  shape,  yielded  absorption  coefficients  less  than  actually  observed 
for  the  10  um  band.  In  addition,  the  absorption  from  wings  of  these 
Lorentzian  line  shapes  in  the  10  um  band  increases  with  increasing 
temperature,  whereas  experiment  has  shown  that  the  actual  continuum 
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1.1  Low  resolution  solar  spectrum  compared  with  laboratory 


spectra  of  atmospheric  gasses  [9]. 


absorption  has  a  negative  temperature  dependence,  that  is,  it  decreases 
with  increasing  temperature  [9], 

This  failure  of  the  Lorentzian  far  wings  to  model  the  continuum  has 
popularized  a  theory  of  continuum  absorption  which  bases  the  phenomenon 
on  water  dimers  [10,11]  or  larger,  more  complex  clusters  or  aggregates 
[12,13]  creating  intense,  unresolved  absorption  bands  throughout  the 
infrared  spectrum.  The  theory  argues  that,  though  the  number  of  these 
complex  aggregates  is  small,  their  effect  could  eas'ly  be  seen  in  the 
window  regions.  Most  importantly,  the  dimer  theory  does  explain  the 
negative  temperature  dependence  of  the  continuum  for  it  is  argued  that 
the  attenuation  is  proportional  to  the  number  of  aggregates  which,  in 
turn,  is  proportional  to  exp(-AH/^T).  where  AH  is  the  change  in  enthalpy 
on  aggregate  formulation  and  is  negative.  Thus,  the  attenuation  will 
show  a  negative  dependence  on  temperature. 

The  inability  of  the  Lorentz  line  shape  to  correctly  predict  the 
magnitude,  temperature,  and  pressure  dependence  of  the  water  vapor 
continuum  has  made  the  dimer  theory  popular.  In  1979,  however,  M.E. 
Thomas  and  R.J.  Nordstrom  at  The  Ohio  State  University  [14]  developed  a 
line  shape  which  is  Lorentzian  out  to  5  cm'l  from  line  center  but  is 
modified  beyond  5  cm*l.  Figure  1.2  shows  the  inability  of  the  Lorentz 
line  shape  to  model  the  water  vapor  pressure  dependence  of  the  nitrogen 
broadened  absorption  coefficient  at  room  temperature  for  the  P(24)  CO? 
laser  line  (940.548  cm'l).  If  only  the  absorption  lines  within  the 
commonly  used  ±  20  cm_l  bound  are  involved  in  the  calculation,  the 
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Figure  1.2  Comparison  of  experimental  and  theoretical  pressure 
dependence  of  water  vapor  continuum  absorption  at 
944.1945  cm-1. 
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Lorentz  shape  falls  far  below  the  actual  data  plot.  But  when  the  bound 
is  extended  to  all  absorption  lines  between  0  and  2300  cm_l,  tne  Lorentz 
line  shape  does  a  much  better  job  of  modeling  the  data.  Both  of  the 
above  calculations  used  the  AFGL  line  listing  [15].  Now  observe  how 
well  the  Thomas-Nordstrom  non-Lorentzian  line  shape  (or  TN  model)  models 
the  pressure  dependence  of  the  absorption  continuum  for  this  10  pm 
window  line. 

Their  study  demonstrated  the  importance  of  far  wing  phenomena  in 
modeling  water  vapor  continuum  absorption  by  first  quantum  mechanically 
deriving  a  model  for  water  vapor  nitrogen  interactions,  and  then  fitting 
the  model  to  data  from  regions  dominated  by  strong  local  line 
absorption  by  fine  tuning  their  line  shape  parameters.  After  the  line 
shape  parameters  had  been  determined  in  this  region,  they  successfully 
applied  the  model  to  the  window  regions  where  continuum  absorption 
dominates.  Data  from  grating  spectrometer  measurements  from  300  cm-1  to 
650  cm*l  [16],  and  HF  [17]  and  CO  [18]  laser  measurements  in  the  3  pm 
and  5  pm  spectral  regions  respectively,  was  used  to  obtain  the  far  wing 
line  shape  parameters  needed  for  the  TN  model.  They  referred  to  these 
parameters  as  the  "A"  and  "G"  parameters. 

After  the  "A"  and  "G"  parameters  were  determined  by  curve  fit  to 
the  above  mentioned  data,  the  line  shape  model  was  applied  to  the  4  and 
10  pm  window  regions  with  great  success  particularly  in  modeling  the 
pressure  dependence  of  the  continuum  absorption  coefficient.  A  negative 
temperature  dependence  was  predicted  by  this  model  but  not  quite  at  the 
rate  actually  observed,  as  Figure  1.3  illustrates.  This  lack  of 


Figure  1.3  Comparison  of  experimental  and  theoretical  pressure 
dependence  of  water  vapor  continuum  absorption  at 
944.1945  cm-1. 
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agreement  is  attributed  to  the  inaccuracies  of  the  interaction 
potentials  and  t.he  perturbation  expansions  of  the  molecular  collision 
Hami Itonian. 

In  addition,  accurate  temperature  dependence  data  does  not  exist 
for  any  of  the  absorption  bands  of  H^O  in  the  infrared  spectral  region. 
With  the  aid  of  such  data,  significant  improvements  and  refinements  in 
the  line  shape  parameters  could  be  made  and  thus  improve  the  model's 
ability  to  predict  the  negative  temperature  dependence  of  the  continuum. 

It  is  the  primary  objective  of  this  study  to  gather  temperature 
dependence  data  in  the  5  urn  spectral  region,  which  is  close  to  the  6.3 
urn  \>2  band  center  and  dominated  by  strong  local  line  absorption,  in 
order  to  lay  the  foundation  for  future  studies  aimed  at  improving  the 
line  shape  parameters  developed  by  Thomas  and  Nordstrom.  Specifically, 
this  data  should  aid  in  improving  the  temperature  dependence  of  these 
far  wing  parameters. 

If  these  far  wing  parameters  are  to  be  used  in  truly  modeling 
continuum  absorption,  they  must  be  independent  of  the  continuum.  For 
this  reason,  data  must  be  taken  in  regions  where  local  line  absorption 
dominates.  Absorption  coefficients  on  the  order  of  1  to  30  km”!  have 
been  obtained  in  the  5.3  urn  region  where  the  continuum  contribution  is 
only  about  0.1  km~l,  from  I  to  2  orders  of  magnitude  smaller  than  the 
local  line  absorbance. 

Data  was  obtained  through  the  use  of  a  multiple  pass  absorption 
cell  and  a  frequency  doubled,  line  selectable  cw  CO2  laser.  The 
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second  harmonic  generation  was  achieved  outside  the  laser  cavity  by  a 
high  efficiency  CdGeAs?  crystal.  The  doubled  CO?  laser  exhibits  a 
spectral  range  from  4.6  to  5.3  um  which  is  similar  to  a  nitrogen  cooled 
CO  laser's  spectral  range.  The  doubled  CO?  laser  system  was  selected 
over  the  CO  laser  because  it  is  much  easier  and  less  expensive  to 
operate  and  more  importantly,  allowed  the  opportunity  to  collect  data  at 
frequencies  different  from  previous  CO  laser  studies,  thus  enhancing  the 
picture  of  the  spectral  line  shapes  in  this  high  absorption  region. 

This  thesis  is  primarily  intended  as  a  progress  report  and  to  aid 
in  future  attempts  to  gather  more  accurate  data  using  the  frequency 
doubled  CO?  laser  probe  system.  The  data  which  will  be  presented  in 
Chapter  IV  needs  to  be  further  refined  and  verified  by  more  data 
acquisition.  This  thesis  offers  recommendations  as  to  how  more  accurate 
data  can  be  obtained  in  the  future. 

Chapter  II  will  qualitatively  develop  the  general  expression  for 
the  absorption  coefficient  given  by  the  Thomas-Nordstrom  model  [14] 
which  was  originally  developed  here  at  The  Ohio  State  University  in 
1976-79.  Then,  Chapter  III  will  discuss  the  experimental  apparatus, 
calibration  tests,  and  procedures  used  to  acquire  the  data  which  will  be 
presented  and  analyzed  in  Chapter  IV.  Both  Chapters  III  and  IV  will 
point  out  the  various  problems  encountered  with  experiment  and  data 
analysis  and  Chapter  V  wi 1 1  offer  a  summary  of  the  project  goals  and 
conclude  with  recommendations  for  future  data  acquisition. 
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CHAPTER  II 


THEORY 


A.  INTRODUCTION 

The  general  expression  for  the  absorption  coefficient  developed  in 
detail  by  Michael  E.  Thomas  and  Robert  J.  Nordstrom  at  The  Ohio  State 
University  [14]  is  qualitatively  developed  in  this  chapter.  Section  B 
develops  the  near  line  center  and  far  wing  line  shapes  separately  and 
Section  C  presents  the  total  line  shape,  the  line  strength  and  the 
final  general  expression  for  the  absorption  coefficient  for  H2O-N2 
gaseous  mixtures. 

Consider  Figure  2.1  which  illustrates  the  propagation  and 
attenuation  of  monochromatic  radiation  by  an  absorbing  medium.  The 
change  in  intensity  dl  per  unit  length  d£  can  be  written  as 

II  =  -ol  (1) 

dt 

where  a  is  the  extinction  coefficient  which  can  further  be  broken  down 
into  two  mechanisms: 

o  *  k  +  s  (?) 
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Figure  2.1  Radiation  transmission  in  a  lossy  medium. 

where  k  is  the  absorption  coefficient  and  s  is  the  scattering 
coefficient. 

This  study  is  restricted  to  infrared  propagation  through  a  gaseous 
medium  (a  <<  X)  thus  making  the  scattering  effects  of  little 
significance.  In  addition,  scattering  is  at  least  a  two  photon  process 
where  absorption  is  a  single  photon  process.  In  the  perturbation 
expansion  of  the  photon  and  absorber  interaction  Hamiltonian,  the  single 
photon  process  will  dominate.  So  the  absorption  mechanism  is  the 
dominant  term  in  the  extinction  coefficient  and  Equation  (?)  can  be 
simplified  to  just  the  absorption  term, 

a  =  k  .  (3) 
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The  absorption  coefficient  indicates  the  probability  of  a  molecular 
interaction  with  an  electromagnetic  field  and  other  molecular  systems 
resulting  in  the  absorption  of  a  photon  and  the  excitation  of  the 
absorbing  molecule. 

Lambert's  law  is  arrived  at  by  the  substitution  of  Equation  (3) 
into  Equation  (1) 

IL.-ki  .  (4) 

dS. 

The  solution  of  this  simple  differential  equation  yields  the 
transmittance  through  the  absorbing  medium  as  Equation  (5)  shows, 

T  =  e'U  (5) 


where 


T 


The  absorption 
attenuation  effects 


coefficient  k(v)  can  be  expressed  as  a  sum  of  the 
of  n  absorption  lines. 


k(v) 


n 

=  ).  Si  ji 


1 


( v-n0;a) 


(6) 


where  S,  is  the  line  strength  of  the  it!l  absorption  line,  and 
j -j  ( v- v0 ;  at)  is  the  line  shape  function  of  the  i^  absorption  line 
centered  at  frequency  v0  and  a  is  the  half  width  at  half  maximum  of 
the  it(l  line. 

The  line  shape  function  is  normalized  by  the  requirement  that 

CO 

J  j(v-v0;a)dv  =  1  (7) 

0 

for  all  i  from  1  to  n.  Thus,  the  absorption  coefficient  can  he  divided 
into  a  line  strength  factor  and  a  line  shape  factor. 

R.  THE  LINE  SHAPE 

1 .  Introduction 

The  two  general  classifications  of  line  shapes  are  homogeneous  and 
inhomogeneous  broadened  line  shape.  Homogeneous  broadening  means  that 
the  line  shapes  have  the  same  basic  shape  for  each  molecule  and  that 
the  shape  observed  for  a  number  of  molecules  will  be  the  same  as  the 
shape  observed  for  one  molecule.  Natural  broadening  and  collision 
broadening  are  examples  of  homogeneous  broadening.  Inhomogeneously 
broadened  lines  are  a  collection  of  frequency  shifted  homogeneously 
broadened  lines  and  the  line  shape  for  a  single  molecule  may  be 
completely  different  than  the  shape  seen  for  a  collection  of  molecules. 
One  example  of  inhomogeneous  broadening  is  Doppler  broadening. 
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Homogeneous  broadening  will  dominate  for  the  pressure  and  temperature 
range  for  which  this  model  is  valid  and  thus  will  be  the  only  type  of 
broadening  considered  in  this  line  shape  development. 

The  homogeneous  line  shape  can  be  broken  up  into  three  regions: 
the  near  line  center  region,  intermediate  region,  and  far  wing  region  as 
shown  in  Figure  2.2  and  Equation  (8), 

=  Wu)  +  jciw<u)  +  Wc^  •  (R) 

The  general  form  of  each  shape  is 

oo 

j  (»)  -  Re  /  dx  exp(-i Au)t)C(t)  (^) 

C  »  o 


Figure  2.2  The  line  shape  function. 
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where  C(t)  is  the  time  correlation  function  for  two  energy  transition 
levels  of  the  absorbing  molecule  perturbed  by  another  molecule  through  a 
collision.  The  lines  labeled  1  and  2  in  Figure  2.3  represent  two 
possible  collision  paths  the  broadening  molecule  might  have  in  relation 
to  the  absorbing  molecule. 

Table  2.1  lists  the  characteristics  of  the  collision  in  each 
region.  The  strong  collisions  typical  of  collision  path  1  greatly  shift 
the  energy  level  of  the  absorbing  molecule  and  thus  are  responsible  for 
the  far  wing  line  shape;  whereas  the  weak  collisions  typical  of 
collision  path  2  only  slightly  perturb  the  absorber  molecular  energy 
level  and  therefore  the  energy  shift  from  line  center  is  very  small  as 
compared  to  the  far  wing  shift. 

The  general  correlation  function  C(t)  is  shown  in  Figure  2.4.  It 
can  be  shown  that  this  function  is  symmetric  about  r  =  0  which  allows 
the  change  in  integration  limits  for  Equation  (9)  resulting  in  the  form 
shown  in  Equation  (ID), 


oo 

j  (<*>)  =  1  Re  /  dr  exp(-iAwi)  exp(iu>  t)C(t) 
c  7n  0 

=  Re  f'!  { g xp ( i  uiq  t)C(  t)  } 


(10) 


Thus,  the  line  shape  function  jc(w)  and  the  function  exp(1u>oT)C(-r)  are 
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TABLE  2.1 


CHARACTERISTICS  OF 
MOLECULAR  COLLISIONS 


R  >  Ra  Very  little  effect  on  absorbing  molecule 


Ra  >  R  >  Rb  Interruption  broadening  [46]  -  near  line 

center  phenomena 

1.  Brief  duration  of  collision 

2.  Time  between  collisions  is  long 
compared  to  duration  of  collisions 

3.  Binary  collisions 

4.  Long  correlation  time 


Rb  >  R  >  Rc  Intermediate  region  -  difficult  to  model 

1.  Many  body  problem 


Rc  >  R  Statistical  broadening  [46]  -  far  winc: 

1.  Binary  collisions 

2.  Duration  of  collision  ~  ^R 

v 

3.  Short  correlation  time 
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Figure  2.3  The  correspondence  between  impact  parameter  and  the  line 
shape  region. 


17 


Figure  2.4  The  correlation  function. 

related  by  the  Fourier  transform.  The  time-bandwidth  product  of  the 
Fourier  transform  theory, 

AtAuj  =1  ( ) 


further  demonstrates  how  the  strong  collisions,  with  short  correlation 
times  at  must  result  in  large  frequency  deviations  from  line  center  Am 
in  order  to  satisfy  Equation  (11). 

The  functional  form  of  the  line  shape  is  arrived  at  by  substituting 
Equation  (10)  into  Equation  (9),  thus  resulting  in  Equation  (12), 


jc(,o)  =  jcFW("}  +  jcIW(w) 

=  Re  F'1  !eja)°T  CFW 

f  Re  E"1  (eju)°T  C 


(t)}  +  Re  F 

Nir(TH 


-1 


|ej-°T 


C[W(T)I 


(1?) 


The  near  line  center  and  far  wing  line  shapes  can  be  determined  with  a 
great  deal  of  work  but  the  intermediate  region  line  shape  is  nearly 
impossible  to  determine.  This  region  will  be  determined  by  a  mixing  of 
the  near  line  center  and  the  far  wing  line  shape  functions  in  Section  3. 

?.  The  Near  Line  Center  Line  Shape  -  Interruption  Broadening 

The  conditions  for  the  near  line  center  phenomena  of  interruption 
broadening  or  impact  broadening  are  listed  in  Table  2.1.  Recause  the 
interaction  time  is  brief  under  these  conditions,  the  correlation 
function  can  be  examined  between  C(t)  and  C(r+dr)  where  dr  <<  t.  That 
is,  the  same  correlation  function  applies  at  t  and  t+dt.  The  result  is 
that  the  correlation  function  appears  as  a  monotonic  decay  during  this 
time  interval  which  can  he  expressed  as  a  simple  exponential  of  the  form 

C( t)  =  exp (-at)  .  (13) 

From  Equation  (13),  the  line  shape  becomes 


JcNLC 


M  =  I 


An 


(14) 


which  is  the  Lorentz  line  shape  and  a  the  halfwidth. 

The  temperature  dependence  of  a  is  difficult  to  derive.  A  second 
order  perturbation  expansion  of  the  molecular  interaction  potentials 
yields  a  temperature  dependence  of  the  form 


a 


«  1 

7T 


(15) 
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However,  the  form  is  not  always  observed  [19,20],  Using  another  method 
which  looks  at  the  molecular  perturbations  strictly  as  phase  modulations 
of  the  eigenstates,  called  the  phase  shift  or  phase  modulation  method, 
leads  to  a  temperature  dependence  of  the  form 

o-i  O6) 

T 

for  the  dipole-dipole  interactions  found  in  self-broadening  and 

a  a _ 1  (17) 

jOT 

for  dipole-quadrapole  interactions  found  in  foreign-broadening.  The 
form  showing  the  pressure  and  temperature  dependence  of  the  halfwidth  is 
shown  in  Equation  (18) 


where 

°o^y^  pb  =  “FOREIGN  BROADENED 

B  =  “SELF  BROADENED 
^FOREIGN  BROADENED 
p  =  absorber  pressure 

=  broadener  pressure 


T  =  temperature 


3.  The  Far  Wing  Line  Shape  -  Statistical  Broadening 

The  far  wing  phenomena,  or  statistical  broadening,  is  characterized 
by  short  correlation  functions  resulting  from  strong  binary  molecular 
interactions.  Due  to  the  strength  of  such  interactions,  no  perturbation 
theory  can  be  applied.  However,  the  phase  shift  method  briefly 
described  in  the  last  section  can  be  used  for  statistical  broadening  to 
determine  the  far  wing  line  shape.  That  is,  by  considering  the  phase 
shift  effects  in  the  eigenstates  to  all  orders,  brought  about  by  the 
strong  molecular  interactions  during  photon  absorption,  the  far  wing 
line  shape  can  be  approximated. 

This  phase  shift  or  phase  modulation  theory  and  another  statistical 
theory  used  by  Fomin  and  Tvorogov  [21]  yield  the  same  leading  term  for 
the  far  wi ng  shape 


'cFW 


(to). 


(A<o)1+3/mj 


(19) 


where  mj  =  6  for  sel f -broadeni ng  and  mj  =  4  for  f orei gn-broadeni ng  by 
N?  or  Oj.  A  more  exact  expression  for  the  far  wing  shape  resulting  from 
the  phase  shift  method  is  shown  in  Equation  (20), 


■’cFW(v) 


=  1 


.6266  \a5  +  -8491  ,^4  j 

|  V-V  I  1  ^  |  V-  V  I '  1  •  7  5 

0  O1 


(20) 
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where 


v  =  wavenumber  at  line  center 
o 

xa6  ■  AIM  q®)1-5  Pa 

xb4  .  A2(v,J)  (M6)1-5  pb 

v  =  vibrational  quantum  number 
J  =  rotational  quantum  number. 

The  frequency  w  (radians/second)  has  been  replaced  by  the  wavenumber 
v  (cm*l).  Xag  and  Xb4  are  very  complicated  functions  of  the  interaction 
potentials  between  H2O  -  H2O  and  H2O  -  N2  interactions,  respectively. 
M(v)  will  be  different  for  every  vibrational  quantum  number  of  the 
absorber  but  is  independent  of  the  rotational  quantum  number  due  to  the 
H2O  -  H?0  interaction  characteristics.  However,  A2(v,J)  does  depend  on 
both  the  vibrational  and  rotational  quantum  numbers.  Since  the 
potential  function  for  the  near  line  center  and  far  wing  phenomena  are 
the  same,  the  j-dependence  of  A2(v,J)  can  be  written  as 


A2(v,J)  =  A2(v)  (19.919)  (^_)1,l?5 


where  0.07  represents  a  mean  value  for  Oq.  Al(v)  and  A2(v)  are  two  of 
the  four  parameters  in  the  TN  model  that  were  determined  experimentally 
because  of  the  complex  nature  of  the  terms  that  make  up  the  parameters 
and  the  uncertainty  of  the  interaction  potentials.  As  stated  in  Chapter 
I,  it  is  the  intent  of  this  study  to  gather  data  that  can  be  used  to 
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improve  upon  the  temperature  modeling  capability  of  the  IN  "A" 
parameters  just  described  and  the  TN  "G"  parameters  which  will  be 
described  shortly. 

Equation  (22)  shows  the  form  that  the  far  wing  shape  takes  on  after 

00 

it  is  normalized  to  one  (i.e.,  /  jcVW(v)dv  =  1). 

o 


JcFW 


. 3 1 98  Xa  5  +  .4334  Xb4 

|Auj1.5+(xa6)3  !Av|1.75+(xb4)7/3 


(22) 


There  is  one  final  correction  that  must  be  made  in  the  far  wing 
line  shape.  The  phase  shift  method  used  to  develop  Equation  (22) 
requires  an  energy  level  shift  or  frequency  modulation  correction  factor 
which  was  determined  by  Fomin  and  Tvorogov  [21].  When  this  level  shift 
correction  is  applied  to  the  far  wing  shape.  Equation  (22)  becomes 


W»>  ■  i 


,  exp(-Ga6(v)/Z;) 


Avj1-5+(Xa6)3 


+  _ exp(-Gb4(  v)/Av)  | 


Av|1.75+(xb4)7/3 


(23) 


where  the  "G"  parameters  are  functions  of  the  vibrational  quantum  number 
of  the  absorber.  The  temperature  dependence  of  the  "G"  parameters  is  of 
the  form  shown  in  Equations  (24)  and  (25), 


(24) 


r 


r,a6(«)  •  Ga(v)  (2«6)-67 
Gb,(»)  ■  Gb(»)  (i«)-67 

These  are  the  remaining  two  of  the  total  four  parameters  which  had  to  he 
experimentally  determined  for  the  TN  model. 

C.  THE  TOTAL  LINE  SHAPE,  THE  LINE  STRENGTH  AND  THE  ABSORPTION 
COEFFICIENT 

The  only  remaining  portion  of  the  line  shape  to  be  determined  is  in 
the  intermediate  region.  By  choosing  band  pass  and  band  stop  filters 
centered  at  v0,  with  slow  fall  off,  a  mixing  of  the  interruption  and 
statistical  line  shapes,  respectively,  can  be  achieved  in  this  region. 
Thus  the  total  line  shape  will  have  the  form 

JCM  =  n^'cNLC(v)P^v^  +  jcFW(v)(1”P(v))  (26) 

where  N  is  the  normalization  constant  determined  by  evaluating  Equation 

(27) , 

oo 

!  d v  j  (v)  =  1  .  (27) 

o 

The  pass  filter  function  was  chosen  to  be  of  the  form  shown  in  Equation 

(28)  and  displayed  in  Figure  2.5, 
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Air  (  cm-' ) 


Figure  2.5  p(u)  versus  Av. 


1  +  .1  cos^.Av  |  Av|  <5 


pH  = 


cm 


Av|  >5  cm" 


(28) 


This  allows  the  interruption  broadening  line  shape  to  dominate  near  line 
center  and  the  band  stop  filter  (l-p(v))  allows  the  statistical 
broadening  line  shape  to  dominate  in  the  far  wings  beyond  5  oH  (50  to 
100  halfwidths  from  line  center). 

One  additioiia1  term  shown  below  must  multiply  the  line  shape 
function  for  v0  <  1000  cm*-. 


v 


(1-e-^c/kT ) 


This  term  accounts  for  the  Maxwel 1 -Boltzman  population  statistics  away 
from  line  center  and  this  term  is  essentially  equal  to  one  for 
u0  >  1000  cm~l.  Therefore,  it  is  not  included  for  that  range  in  order 
to  improve  the  absorption  coefficient  computation  time. 

The  line  strength  is  derived  by  considering  the  transition  rates 
from  both  the  lower  to  upper  states  and  the  upper  to  lower  states  and 
can  be  shown  to  have  the  following  form. 


S 


£u 


Va  M 

TPa  T 
o 


E 

e 


(T-T0) 
1  kTT0 


(1-e 

(1-e 


-vQhc/kT 
-v  hc/kT 


L 

°) 


(29) 


where 


Stu  =  the  1ine  strength  at  standard  pressure  and  temperature 

T0  =  296°k 

P  =1  atm 
ao 

The  total  absorption  coefficient  for  H2O  -  N2  mixtures  is  given  in 
Equation  (30), 

k(v)  =  l  k-j  ( v)  =  l  Stu(v)  jc(v)  .  (30) 

i  iu 

The  summation  is  over  all  water  vapor  transitions.  Using  the  form  shown 
in  Equation  (26),  the  HpO  -  N2  absorption  coefficient  becomes 
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(31) 


k(v)  =  I  Siu  N  [  jCN|_c(  V)P  ( v)  +  JcFw(  v)  ( 1-P  ( v) )  1 
2u 

where 


ToPa 

*u  TpT 
o 


E  (T-To) 
e  1  kTT0 


(i-g 

(l-e 


- v  hc/kT 
0 

-v  hc/kT 
o 


L 

°) 


jcNCL^  = 


i  (1-e-^c/kT) 

*  vo  - v0hc/kT 
(l-e 


)((Av)2+(aP)2) 


1  of 

*  (Au)2+(ar)2 


v  <  inoo  cm'1 
vQ  >  1000  cm'1 


1  (l-e-^cAT) 

»  ,^-v.cAT, 


.3198^6  e'Ga6’'“”  _4334Jm  e-Gb4/«» 


|  Av| 1  -5+x3 

3b 


|  av| 1.75+x7^3 
-1 


jcFW^' 


vq  <  1000  cm 


•jiqqi  .  -Gs4/av 

■3198Xag  e  .4334X^4  e 


|  Avj 1  *^+^a6 


AvI  1.75+x7/3 
b4 


vq  >  1000  cm" 


The  normal ization  constant,  N,  is  determined  by 


N_1  =  I  JcNLC ( v)  P(v)dv  +  /  jcFW(v)  (l-p(v))dv 
0  0 


(32) 
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CHAPTER  III 


EXPERIMENTAL  APPARATUS,  TESTS  AND  PROCEDURES 

A.  INTRODUCTION 

This  chapter  describes  the  apparatus  used  for  this  study,  the 
calibration  tests  performed  on  the  data  acquisition  system  and  the 
procedures  used  in  sample  handling  and  data  acquisition.  The  apparatus 
description  begins  with  the  CO2  laser  followed  by  a  description  of  the 
second  harmonic  generation  set  up,  the  optics  external  to  the  absorption 
cell,  the  absorption  cell  and  the  data  acquisition  system. 

B.  THE  EXPERIMENTAL  APPARATUS 
1.  The  CO2  Laser 

The  laser  used  in  these  experiments  is  an  axial  flow  cw  CO2  laser 
which  is  grating  tunable  for  single  line  operation.  The  output  power 
for  the  laser  lines  used  in  the  study  varied  from  2.5  to  3.0  watts.  The 
grating  can  be  tuned  over  both  the  10.6  urn  and  9.6um  bands  with  output 
power  greater  than  1  watt  for  25  and  22  lines,  respectively.  Table  3.1 
shows  the  typical  observed  power  for  each  of  the  lines. 
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TABLE  3.1 


TYPICAL  OUTPUT  POWER  FOR  EACH  CO2  ‘.LASER  LINE 
00°1  -  02°0  BAND  9.6  um 


LINE 

POWER  (WATTS) 

LINE 

POWER  (WATTS) 

R(26) 

1.1 

P(10) 

1.8 

R(24) 

1.9 

P  ( 12 ) 

3.2 

R(22) 

2.8 

P  ( 14 ) 

3.1 

R  (20) 

2.6 

P(16) 

3.3 

R(18) 

2.8 

P(18) 

3.4 

R(16) 

2.1 

P(20) 

4.2 

R(14) 

2.3 

P(22) 

4.0 

R(12) 

1.2 

P  ( 24 ) 

4.0 

R(10) 

1.3 

P  ( 26 ) 

3.2 

R(8) 

0.9 

P(28) 

2.9 

P  (30) 

3.0 

P  (32) 

1.5 

P  ( 34 ) 

0.9 

P  (36) 

1.1 

00°1  - 

10°0  BAND  10.6  um 

LINE 

POWER  (WATTS) 

LINE 

POWER  (WATTS) 

R  (32) 

0.9 

8(6) 

0.7 

R  (30) 

1.5 

8(8) 

1.2 

R  (28) 

2.2 

8(10) 

1.9 

R(26) 

3.2 

8(12) 

2.1 

R  (24) 

2.9 

8(14) 

3.1 

R  (22) 

3.1 

8(16) 

3.0 

R  (20) 

3.8 

8(18) 

3.0 

8(18) 

3.0 

P(20) 

3.8 

8(16) 

2.5 

P(22) 

4.6 

8(14) 

3.0 

8(24) 

3.0 

8(12) 

2.8 

8(26) 

3.7 

8(10) 

2.0 

P(28) 

2.5 

8(8) 

1.1 

8(30) 

2.0 

8(32) 

1.1 

8(34) 

0.9 
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The  resonant  cavity  for  the  laser  that  is  shown  schematically  in 
Figure  3.1  consists  of  a  flat  grating,  an  adjustable  iris,  a  ZnSe 
Brewster  angle  window  and  a  ZnSe  output  window  which  is  60%  reflecting 
at  10.6  ym  and  has  a  5  meter  radius  of  curvature.  The  output  mirror  is 
mounted  on  an  adjustable  piezoelectric-dri ven  mount  which  provides 
sufficient  length  tuning  of  the  cavity  for  optimal  power  output  and  open 
loop  stability. 

The  resonant  cavity  stability  condition  is 

0  <  (1  -  2L)  (1  -  2L)  <  l  (33) 

where  L  =  1.5  meters,  the  cavity  length;  =  °°,  the  radius  of  curvature 
of  the  flat  grating,  and  R2  =  5  meters,  which  is  the  radius  of  curvature 
of  the  output  mirror.  Indeed,  this  cavity  meets  this  condition  as  shown 
in  Equation  (34), 


0  <  0.40  <  1  .  (34) 

Calculations  of  the  beam  waist  at  the  grating  and  the  spot  size  at  the 
output  mirror  are  shown  in  Equations  (35)  and  (36), 
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where 


9  "  1  '  i?2  '  0-7 

L  =  1.5  m 

X  =  10.6  pm 

The  iris  was  required  in  the  cavity  to  insure  single  mode 
operation.  The  grating  is  tuned  with  a  micrometer  adjustment  and  has  a 
relatively  large  dispersion  which  permits  precise  and  reliable  selection 
of  the  desired  line. 

The  laser  is  an  axial  flow  type  of  system  containing  a  gas  mixture 
of  He-N2*C02  with  a  partial  pressure  ratio  of  80:10:10,  respectively. 

The  gas  pressure  at  the  inlet  is  typically  11.5  Torr.  Tap  water  flows 
through  an  external  jacket  and  serves  to  cool  the  gas  flow.  The  laser 
discharge  is  excited  by  a  Universal  Voltronics  model  BAC  16-50  DC  power 
supply  connected  in  series  to  a  forced  air  cooled  700  Kq  ballast 
resistor/capacitor  bank.  Typically,  the  DC  voltage  measured  at  the 
source  is  13  KV  with  a  9  mA  current. 

The  laser  tube  is  supported  by  a  1.7  meter  long  Invar  structure 
mounted  to  the  steel  base  plate  at  the  grating  end  only,  with  the  other 
end  resting  in  a  sliding  mount  which  allows  for  thermal  linear  expansion 
of  the  structure. 
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2.  Second  Harmonic  Generation 

The  second  harmonic  generation  was  achieved  by  mounting  a  CdGeAs2 
crystal  on  the  cold  finger  of  a  Jam’s  research  dewar  model  DT,  cooling 
it  down  to  77K  with  liquid  nitrogen,  and  then  pumping  th'  crystal  with 
approximately  3  watts  of  cw  radiation  from  the  previously  described  COj 
laser.  Figures  3.2  and  3.3  show  the  crystal  dewar  and  the  r.C>2  laser, 
respectively. 

The  CdGeAs2  crystal,  which  is  3  x  5.5  mm  in  cross  section  and  6.4 
mm  long,  was  loaned  to  us  by  Dr.  Dennis  Killinger  and  Dr.  Norman  Menyuk 
of  The  Massachusetts  Institute  of  Technology  Lincoln  Laboratory.  This 
crystal  is  one  of  the  most  useful  second  harmonic  generation  materials 
for  the  infrared  spectral  region,  since  it  has  one  of  the  highest 
non-linear  coefficients  known  and  a  useful  transmission  range  from  3  urn 
to  17  urn. 

The  optical  absorption  is  decreased  by  cooling  the  crystal,  with 
Lincoln  Laboratory  reporting  0.1  and  0.4  cm"!  at  10  and  5  urn, 
respectively,  having  been  measured  for  their  best  samples  [22].  Our 
crystal  shows  considerably  higher  absorption  but  generates  the  desired 
amount,  of  second  harmonic  energy. 

It.  was  cut  for  Type-I  phase  matching,  with  a  propagation  angle  of 
32.5°  from  the  optical  axis.  Figure  3.4  shows  how  the  fundamental  must 
be  polarized  and  propagated  into  the  crystal  for  phase  matching  with  the 
desired  second  harmonic.  The  second  harmonic  emerges  from  the  crystal 
polarized  in  the  ordinary  direction,  that  is,  the  polarization  is 
perpendicular  to  the  optical  axis,  whereas  the  incoming  fundamental  is 
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Figure  3.3  1.5  meter  C02  laser. 


3.4  Polarization  of  input  and  output  fields  for  second  harmonic 
Type-I  phase  matching. 


1 


polarized  in  the  extraordinary  direction.  One  of  the  important 
properties  of  all  nonlinear  crystals  is  that  their  index  of  refraction 
varies  as  a  function  of  both  frequency  and  direction  of  propagation  of 
the  polarized  radiation.  These  properties  of  birefringence  and 
dispersion,  respectively,  are  fundamental  to  second  harmonic 
generation  [231. 

Figure  3.5  illustrates  the  elliptical  characteristics  of  the  index 
of  refraction  as  a  function  of  the  spatial  orientation  of  the  direction 
of  polarization  for  an  electromagnetic  wave  propagating  through  a 
typical  nonlinear  crystal.  For  the  extraordi nary  direction  of 
po1 jri zat i on ,  an  e-ray,  the  index  changes  elliptically  between  the  value 
of  the  ordinary  index  n0,  when  the  wave  normal  is  parallel  to  the 
optical  axis,  and  the  ext raordi nary  index  ne,  when  the  wave  normal  is 
perpendicular  to  the  optical  axis.  Phase  matching  and  therefore  second 
harmonic  generation  is  achieved  when 

n*  =ne(h)  (37) 

2w  w 

where  n^w  =  second  harmonic  o-ray  index,  and 


ne  -  fundamental  e-ray  index  as  a  function  of  9. 

w 


Because  the  index  of  the  fundamental  e-ray  can  be  varied  by  changing  the 
angle  between  the  optical  axis  and  the  wave  normal,  we  are  able  to 
propagate  the  wave  at  an  angle  9  to  the  optical  axis  such  that  the 
refractive  index  for  the  input  is  exactly  equal  to  the  refractive  index 
for  the  second  harmonic.  This  angle  is  approximately  9  =  35°  for  10.6 
urn  and  thus  is  the  reason  for  the  crystal  being  cut  as  shown  in 
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Figure  3.5  Illustration  of  the  elliptical  characteristics  of  the  index 
of  refraction  as  a  function  of  the  spatial  orientation  of 
polarization  for  an  electro-magnetic  wave  propagating 
through  a  typical  nonlinear  crystal. 
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Figure  3.2.  With  the  crystal  cut  in  this  way,  the  phase  matched 
fundamental  wave  normal  is  nearly  perpendicular  to  the  crystal  face. 


The  second  concern  in  generating  second  harmonics  is  the 
fundamental  beam  intensity  inside  the  crystal.  It  can  be  shown  [23] 
that  there  is  an  optimum  beam  waist  size  for  a  given  wavelength  and 
crystal  length  in  order  to  achieve  the  maximum  second  harmonic  output 
power.  This  optimum  beam  waist  size  is  obtained  by  setting  L  =  6.4  ran 
and  X  =  10.6  urn  in  Equation  (38), 


2tt 


62.14  uro 


(38) 


The  beam  was  focused  in  the  crystal  with  a  30  cm  radius  of  curvature 
spherical  mirror.  The  spot  size  on  the  mirror  was  calculated  to  be 
ui  =  4.15  urn  so  the  beam  waist  at  the  focus  was  calculated  to  be 
oj0  =  121.95  m  using  Equation  (39) 

or  =  11  =  121.95  ym  (39) 

0  mu) 


where 

f  =  R/2 
R  =  30  cm 

This  is  roughly  twice  the  optimum  beam  waist;  however,  the  second 
harmonic  power  obtained  with  this  intensity  was  sufficiently  high 
for  our  purposes.  In  addition,  focusing  with  the  larger  beam 
waist  lowers  the  possibility  of  thermally  damaging  the  crystal. 
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3 .  The  c xtereal  Optical  System 

The  entire  optical  system,  external  to  the  absorption  cell,  is 
shown  in  Figure  3.6.  Most  of  the  10.6  urn  radiation  is  attenuated  by  a 
L i F2  window  mounted  to  the  output  window  of  the  crystal  dewar.  This 
filtering,  together  with  the  sapphire  windows  on  the  detectors  and  their 
frequency  response  to  10.6  urn  insure  that  the  radiation  measured  during 
these  experiments  was  the  5.3  um  second  harmonic. 

Other  components  of  the  external  optical  system  include  a  Hughes 
Helium-Neon  laser  model  322H-PC  (HeNe  #1  of  Figure  3.6),  which  was  bore 
sighted  with  the  CO2  laser  and  allowed  precise  alignment  of  the  CO2 
laser  with  the  crystal  and  absorption  cell  optics,  and  a  Spectra  Physics 
Helium-Neon  laser  model  125  (HeNe  §2)  which  served  as  the  primary 
alignment  laser  for  the  absorption  cell  optics  and  detectors.  The  two 
irises  were  used  to  colinearize  the  two  alignment  lasers  and  thus  align 
the  second  harmonics  with  the  absorption  cell  and  the  detectors. 

4.  The  10.785  Meter  Multipass  Absorption  Cell 

The  24  inch  diameter  by  10.785  meter  long  multipass  absorption  cell 
used  for  this  experiment  was  designed  by  Professor  Edward  K.  Damon  [15] 
and  the  unique  optics  system,  shown  in  the  top  portion  of  Figure  3.6, 
which  couples  the  laser  and  the  interferometer  to  the  optics  of  the 
absorption  cell  was  designed  by  Dr.  Robert  J.  Nordstrom.  The  10.875 
meter  cell  design  is  a  significant  improvement  over  the  old  15.24  meter 
cell  previously  described  [25,26]  and  used  at  The  Ohio  State  University, 
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as  wel 1  as  the  standard  White  type  cell  design  [27].  The  improved 
design  allows  for  longer  path  lengths  with  improved  spatial  stability 
and  reduced  aberration  effects. 

The  internal  mirrors  have  a  relectivity  of  98.67,  at  10.6  pm  and  a 
radius  of  curvature  of  10.785  meters.  Figure  3.7  shows  the  spot 
arrangement  on  the  field  mirror  of  the  cell.  Spot  5  would  typically  be 
the  output  spot  on  a  White  type  system  but  on  this  cell  it  is  fed  back 
into  the  optical  system  generating  two  additional  rows  of  spots.  Path 
lengths  up  to  1.8  Km  are  possible  with  laser  sources  and  lengths  of  1.0 
Km  have  been  achieved  with  broadband  sources.  The  mirrors  are  adjusted 
by  microprocessor  controlled  stepper  motors  thus  allowing  the  path 
lengths  to  be  easily  changed. 

Calculation  of  the  path  length  is  simply  determined  by  considering 
that  each  spot,  including  spot  5  and  the  output  spot,  represents  a  round 
trip  length  of  21.570  meters.  An  additional  1.5  meters  is  added  to 
account  for  the  path  lengths  of  the  input  and  output  periscope  optics. 
The  path  length  used  for  this  study  was  130.9  meters,  where  a  column  of 
four  spots  was  seen  on  the  field  mirror  plus  spot  5  and  the  output  spot. 

Other  features  of  the  cell  include  a  mounting  system  that  rigidly 
couples  the  external  optics  to  the  cell  but  isolates  both  the  external 
optics  and  the  cell  from  any  building  vibration.  The  cell  walls  are 
constructed  from  304  stainless  steel  with  a  very  smooth  interior  to  keep 
absorbance  effects  to  a  minimum.  Buffer  gases  are  introduced 
uniformly  along  the  cell  allowing  for  much  faster  sample  mix  times 
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than  those  of  the  old  15.24  meter  cell.  The  vacuum  system  is  a  Kinney 
model  KTC-112  compound  vacuum  pump  which  is  capable  of  evacuating  the 
entire  cell  in  about  5  minutes.  Another  very  important  feature  of  this 
cell  is  that  it  can  be  temperature  controlled  from  -60°C  to  +60°C  with  a 
uniformity  of  2  C  over  the  entire  cell  length. 

5.  The  Computer  Controlled  Data  Acquisition  System 

The  data  acquisition  system  used  for  this  experiment  consists  of  an 
IMSAI  Model  8080  mi crocomputer ,  a  teletype,  a  high  speed  24  channel 
analog  to  digital  converter,  two  DC  amplifiers,  two  lock-in  amplifiers, 
a  pre-amp  and  filter  unit,  and  two  InSb  detectors. 

The  two  channels  of  the  A/D  are  sampled  1024  times  each  with  the 
channels  being  sampled  10  msec  apart  with  a  50  msec  delay  between  each 
pair  of  readings.  The  1024  samples  are  divided  into  four  groups  of  256 
each  and  after  a  group  of  256  samples  has  been  taken  a  printout  is 
made  of  the  average  and  standard  deviation  for  each  channel  as  well  as 
the  ratio  of  the  two  channels,  i.e.,  Si gnal /Reference.  This  grouping 
allows  the  experimenter  to  observe  any  trends  in  the  data.  The 
averaging  tends  to  filter  out  any  signal  variations  due  to  120  Hertz 
line  ripple,  room  vibration,  or  laser  power  fluctuations. 

Figure  3.8  shows  a  block  diagram  of  the  data  acquisition  system. 

The  A/0  unit  is  a  12  bit  Altair  A/D  with  2.5  mV  resolution,  measuring 
from  0  to  9.9974  volts.  Fctron  model  R516-6A  DC  amplifiers  were  used  to 
amplify  the  analog  signals  coming  out  of  the  lock-in  amplifiers. 
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Their  output  signal  ranges  from  0  to  1  volts  so  the  DC  amps  were  set  to 
a  gain  of  10  to  boost  the  signal  levels  to  between  0  and  10  volts  to  the 
A/D  unit.  The  calibration  of  the  system  is  discussed  in  Section  C. 


Princeton  Applied  Research  model  128A  lock  in  amplifiers  were  used 
in  this  data  acquisition  system  to  synchronously  detect  and  integrate 
the  signal  from  each  detector.  A  radiation  chopper.  Laser  Precision 
model  CTX-534,  provided  100%  modulation  of  the  second  harmonic  into  the 
detectors  and  has  its  external  synchronizing  line  connected  to  the  phase 
reference  input  of  the  lock  in  ampifiers. 

Specifically,  the  lock  in  amplifiers  provide  accurate  voltage 
measurement,  for  signals  corrupted  by  broad  hand  noise  and  power  line 
pick  up,  for  example.  An  extremely  narrow  band  synchronous  detector, 
whose  passband  center  is  locked  in  with  the  chopper  frequency  provides 
an  excellent  signal  to  noise  ratio.  The  input  selectivity  range  can  be 
varied  from  1  V  to  250  mV  full  scale.  The  analog  output  has  an 
adjustable  low  pass  f i Iter  which  can  be  used  to  damp  out  fluctuations 
from  the  synchronous  detector  stage. 

1  nSti  detectors  are  excellent  for  detection  of  small  signal 
radiation  in  this  spectral  region  for  their  frequency  response  peaks  at 
about  4.6  um  which  is  close  to  the  5.3  urn  second  harmonic  radiation. 

The  signal  detector,  a  Judson  model  J-10  nitrogen  cooled  InSb  detector, 
is  photovoltaic  and  thus  requires  no  biasing  circuitry.  Its  output  was 
connected  directly  to  the  signal  channel  lock  in  amplifier.  The 
reference  detector,  an  Infrared  Associates  model  DNSL-  18  Hrs.  nitrogen 


rooloH  inc.h  dptect.or,  is  also  photovoltaic.  However,  because  of  the 
very  low  intensity  of  the  reference  beam  after  reflection  from  a  BaF2 
wedge,  it  was  found  that  a  better  signal  to  noise  ratio  could  be 
obtained  by  connecting  the  detector  to  a  Nicolet  FT1R  pre-amp  and  filter 
unit  and  then  to  the  reference  channel  lock-in  amplifier. 

C.  TESTS 

1.  DC  Amplifier  and  A/D  Calibration 

The  Ectron  DC  amplifiers  and  the  Altai r  A/D  were  tested  and 
calibrated  by  feeding  the  output  from  an  Analogic  model  AN  3100  voltage 
standard  with  a  resolution  of  0.1  millivolts  into  each  DC  amplifier  and 
incrementing  the  voltage  from  0.05  to  0.95  volts  by  0.05  volt  steps. 

The  output  from  the  amplifiers  was  read  at  each  step  by  the  A/D  and 
computer.  Since  the  gain  of  the  amplifiers  was  10  the  theoretical 
output  was  from  0.5  volts  to  9.5  by  0.5  volt  increments.  This 
theoretical  value  was  compared  with  the  value  actually  read  and  the 
error  (=output  -  theoretical)  was  plotted  and  is  displayed  in  Figure 
3.9.  This  error  is  very  small  and  averages  to  less  than  0.4%  across  the 
entire  voltage  range.  Since  the  signal  channel  reading  is  divided  by 
the  reference  channel  and  the  signal  and  reference  calibration  curves 
track  with  each  other  so  closely,  this  error  is  almost  completely 
cancelled  out  in  the  data  measurements. 


ERROR  (mV) 
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CALIBRATION  OF  DATA  SYSTEM 


Figure  3.9  Calibration  data  for  the  signal  and  reference  channels 
The  error  in  mV  is  plotted  versus  theoretical  output. 
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2.  The  Lock-in  Amplifiers 

Both  of  the  lock-in  amplifiers  had  been  recently  calibrated  and 
tested  in  a  previous  study  [26]  and  so  only  a  few  non-rigorous  tests 
were  conducted.  A  continuity  check  was  made  between  the  different 
sensitivity  scales,  a  rough  calibration  check  vas  made  using  a 
sinusoidal  input,  and  the  lock-in  amplifiers  were  also  checked  for 
reference  detector  phase  drift.  All  of  the  above  mentioned  tests 
indicated  that  the  amplifiers  were  still  performing  quite  well. 

D.  EXPERIMENTAL  PROCEDURES 

1.  Introduction 

This  section  discusses  the  general  procedures  used  to  align  the 
optical  system,  get  the  proper  amount  of  sample  in  the  absorption  cell, 
and  to  acquire  the  data.  In  addition,  various  technical  difficulties 
encountered  throughout  the  experiment  will  be  discussed. 

2.  Optical  Alignment 

Much  of  this  discussion  refers  to  Figure  3.6.  Initially  the  large 
HeNe  #2  shown  in  Figure  3.6  was  used  to  align  the  absorption  cell  optics 
and  to  set  the  number  of  passes  through  the  cell  to  six  for  a  path 
length  of  130.9  meters.  Then  HeNe  #2  was  used  to  align  the  signal  and 
reference  detectors  for  optimum  response.  Next,  with  the  CdGeAs2 
crystal  temporarily  moved  out  of  the  optical  path,  the  HeNe  #1  laser 
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(Figure  3.6)  which  is  boresighted  with  the  CO2  laser  was  colinearized 
with  HeNe  #2  through  the  use  of  Iris  1  and  2  and  adjustment  of  mirrors 
A  and  8.  At  this  point  the  detector  alignment  was  rechecked  using 
HeNe  #1. 

The  crystal  was  then  moved  back  into  the  optical  path  using  HeNe  #1 
to  target  the  crystal.  The  crystal  had  to  be  targeted  with  liquid 
nitrogen  already  in  the  crystal  dewar  because  the  rod  that  the  crystal 
was  mounted  on  inside  the  dewar  would  contract  by  approximately  2 
millimeters  during  the  initial  cooling  down  to  the  operating  temperature 
of  77°K.  Once  the  crystal  was  in  place,  the  CO2  laser  was  turned  on  and 
the  second  harmonics  were  tuned  in  by  translating  and  rotating  the 
entire  crystal  dewar  assembly.  In  addition,  the  crystal  could  be 
vertically  adjusted  by  a  pedestal  which  the  dewar  rested  on.  This 
pedestal  had  to  be  abandoned  later  due  to  the  fact  that  thermal 
expansion  and  contraction  of  the  dewar's  outer  skin  vertically 
disoriented  the  crystal  and  so  the  dewar  had  to  be  suspended  to  null 
these  effects.  After  a  sufficient  amount  of  second  harmonics  had  been 
tuned  in  by  crystal  reorientation,  mirrors  A  and  B  were  again  used  to 
recheck  and  correct  the  detector  alignment  for  optimum  response. 

3.  Background  Measurements 

After  establishing  a  sufficient  amount  of  second  harmonic  (i.e., 

2.5  mv  reading  from  the  signal  detector),  a  100%  transmittance  level  or 
background  would  be  measured.  Equation  (39)  shows  the  formula  used  to 
obtain  the  transmittance  measurements. 
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(40) 


y  -  1 S  (  n  )  / 1  R  S 

T  =  transmittance 

Is(n)  =  attenuated  intensity  level  due  to  the  sample 
Ir$  =  reference  detector  signal  level  during  sample 
measurement 

Ig(n)  =  unattenuated  intensity  level  (signal  background) 

I 3R  =  reference  detector  signal  level  during  signal  background 
measurement 

n  =  number  of  traversals. 

Thus,  as  Equation  (40)  indicates,  a  Ig(n)/lgg  background  must  be  known 
in  order  to  obtain  the  transmittance. 

A  stable  background  is  very  important  for  obtaining  accurate  data 
but  unf ortunately  the  background  measurements  for  this  experiment  were 
quite  often  unstable.  The  instability  was  tolerable  during  data 
acquisition  at  the  first  two  temperatures  of  295°K  and  310°K  but  became 
too  unreasonable  to  work  with  during  data  acquisition  at  the  third 
temperature  of  323°K.  Thus,  the  data  gathered  for  the  first 
temperatures  can  be  taken  with  a  fair  amount  of  confidence  but  the  data 
at  T  =  323°K  is  in  need  of  much  verification. 

The  source  of  the  instability  was  never  fully  determined  although  a 
couple  of  revealing  observations  were  made  near  the  end  of  the  study 
when  the  background  could  not  be  stabilized  at  all.  First  it  was 
discovered  that  a  stable  background  could  not  even  be  obtained  by  just 
using  the  HeNe  #1  laser  as  the  source  (independent  of  the  COg  laser  and 
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the  crystal).  When  the  Judson  J - 10  signal  detector  was  replaced  by  a 
HgCdTe  (12  urn  peak)  detector,  a  stable  background  was  then  able  to  be 
measured  using  HeNe  #1.  However,  when  just  the  CO2  laser  was 
reintroduced  into  the  system  (i.e.,  no  second  harmonic  generation)  with 
the  new  signal  detector,  the  background  was  again  unstable  although  not 
nearly  as  much.  This  seemed  to  indicate  that  the  CO2  laser  was  a  source 
of  instability  also;  although  it  was  never  determined  how  it  could  cause 
such  instability.  In  any  case,  a  sufficient  amount  of  data  was  obtained 
at  the  two  lower  temperatures  when  the  system  was  relatively  stable. 

4.  Sample  Handling 

The  water  samples  were  loaded  into  the  cell  by  first  filling  a 
pyrex  glass  container  which  was  vacuum  coupled  to  the  side  port  of  the 
absorption  cell,  and  then  heating  the  container  with  a  heat  gun,  thereby 
allowing  the  water  to  evaporate  into  the  cell.  The  desired  partial 
pressure  of  the  water  was  established  by  measuring  into  the  sample 
container  the  volume  of  water  required  to  get  the  desired  pressure  for 
some  known  temperature  and  volume  of  the  cell  in  accordance  with  the 
ideal  gas  law.  Table  3.2  lists  the  water  volumes  used  for  the  various 
sample  pressures  and  temperatures  investigated.  Once  the  sample  was  in 
the  cell,  it  was  nitrogen  broadened  to  a  total  pressure  of  740  Torr. 

This  pressure  was  chosen  so  that  the  cell  would  have  a  slight  negative 
pressure,  for  which  the  cell  seals  were  designed,  and  thus  a  stable 
total  pressure  could  be  maintained  throughout  the  transmittance 
measurement. 
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FABLE  3.2 


H20  SAMPLE  VOLUMES  REQUIRED  TO  ESTABLISH  PARTIAL  PRESSURE 
OF  5,  10,  AND  15  TORR  AT  TEMPERATURES  OF 
295°K,  310°K  AND  323°K 


■ 

H20  Partial 
Pressure  (Torr) 

T  =  295°K 
Volume  (m2) 

T  =  310°K 
Volume  (m2.) 

T  =  323°K 
Volume  (mt) 

5 

16.8 

16.0 

15.3 

10 

33.6 

31.9 

30.6 

15 

50.3 

47.9 

45.9 

5 .  Transmittance  Measurements 

Once  the  nitrogen  broadened  sample  was  in  the  cell  and  the  desired 
temperature  established,  several  readings  were  taken  for  each  of  the 
three  laser  lines  of  interest.  If  the  readings  appeared  to  be  slightly 
unstable,  measurements  were  taken  until  a  stable  plateau  was  observed 
and  then  the  cell  quickly  pumped  down  to  obtain  a  background 
measurement.  On  the  other  hand,  if  the  readings  appeared  fairly  stable 
for  a  10  Torr  sample  partial  pressure,  for  example,  the  cell  would  be 
pumped  down  to  just  370  Torr  and  refilled  with  nitrogen  to  740  Torr. 
This  would  cut  the  partial  pressure  of  the  sample  to  5  Torr  and  another 
set  of  transmittance  readings  would  be  obtained.  In  any  case,  a  final 
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pmpt.y  cell  background  was  always  measured.  The  absorption  coefficients 
were  then  determined  using  Equation  (40)  and  Equation  (41)  where  two 
absorption  coefficient  measurements  would  be  determined. 


K  =  tn  T  (km'1) 

7TH5S 

where 

T  =  transmittance 
.1309  =  path  length  in  km. 


(41) 


The  first  measurement  used  the  transmittance  based  on  the  first 
background  obtained  for  introduction  of  the  sample  into  the  cell  and  the 
second  used  the  transmittance  based  on  the  final  background 
reading.  These  two  readings  make  up  the  end  points  of  the  error  bars 
shown  in  the  data  plots  (Figures  4.1  through  4.9).  Although  the 
coefficient  obtained  by  using  the  final  background  reading  may  in  some 
cases  be  more  accurate  due  to  the  shorter  time  interval  between  sample 
and  background  readings,  the  same  weighting  was  given  to  each  end  point 
when  the  curve  of  best  fit  was  determined. 
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CHAPTER  IV 


DATA  ANALYSIS 


A.  INTRODUCTION 

Measurements  of  nitrogen-broadened  water  vapor  absorption  made  at 
temperatures  of  296°K,  310°K  and  323°K  are  presented  in  this  chapter  for 
three  frequency  doubled  CO2  laser  lines.  The  first  section  discusses 
the  line  selection  criteria  used  in  selecting  the  investigated 
frequencies  as  well  as  a  word  about  the  temperature  selection.  Then  the 
data  is  presented  and  compared  to  the  theoretical  absorption  predicted 
by  the  TN  model. 

B.  SELECTION  OF  THE  INVESTIGATED  FREQUENCIES,  TEMPERATURES 

AND  PRESSURES 

This  study  investigated  the  nitrogen  broadened  water  vapor 
absorption  for  three  COj  laser  lines: 

1)  2 v[R ( 12 ) ]  =  1941.0958  cm'1 2 3 

2)  2v[R (14)]  =  1943.8618  crW 

3)  2v[R (18)]  =  1949.2452  cm’1  . 
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The  above  lines  are  from  the  (00°1  -10°0)  10.6  urn  band  of  the  CO2  laser 
and  were  frequency  doubled  by  the  CdGeAs2  crystal  described  in  Chapter 
III. 

Because  the  primary  intent  of  the  study  was  to  observe  the  line 
shape  behavior  of  local  lines  where  the  continuum  absorption  was  of 
little  significance,  these  lines  were  chosen  on  the  basis  of  their 
locations  relative  to  the  line  centers  of  local  absorption  lines.  Due 
to  a  lack  of  740  torr  nitrogen  broadened  spectral  data  in  this  region, 
the  selection  was  based,  in  part,  on  the  pure  water  vapor  spectra  shown 
in  Figures  4.1a  and  4.1b.  Table  4.1  indicates  the  different 
temperatures  and  pressures  for  which  each  of  the  four  rows  of  spectra 
were  recorded.  The  second  row  (Row  B)  best  approximates  the  pressure 
and  temperature  range  of  the  study.  Observe  that  the  2v[R(12)]  and  the 
2 v[R ( 14 ) ]  lines  lie  somewhere  in  the  intermediate  to  near  line  center 
regions  of  the  local  absorption  lines.  The  intermediate  wing  of  the 
absorption  line  shape  has  the  least  known  about  it  theoretically,  as  was 
mentioned  in  Chapter  II,  and  little  is  known  about  this  region 
experimentally  as  a  function  of  temperature.  Therefore,  these  two  laser 
lines  are  of  great  interest. 

Pressure  broadened  spectra  in  this  region  does  need  to  be  recorded 
to  verify  the  position  of  these  two  lines  relative  to  the  neighboring 
local  absorption  lines  because  it  is  uncertian  just  exactly  where  these 
two  laser  lines  sit  on  the  local  line  shapes.  The  2v[R(12)l  = 

1941.0958  cm~l  line  is  0.5332  cm*l  away  from  a  strong  line  at 
1941.6290  cm‘1.  If  we  assume  that  the  pressure  broadened  line  width 


TABLE  4.1 


EXPERIMENTAL  CONDITIONS  FOR  THE  PURE  H20  SPECTRA 
SHOWN  IN  FIGURES  4.1a  AND  4.1b 


(PATH  LENGTH  =  821  METERS) 


ROW  OF  T  H20  PRESSURE 

FIGURE  (°k )  (TORR) 


A 

284.4 

5 

B 

295.3 

10 

C 

307.1 

20 

D 

320.2 

40 
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2  v  r  l»  ;  1 2  >  ]  2  c  [  R  ( 1 4  )  ]  2  1C  [  R  :  16  )  j  2  ic  [  R  ( 1 8 1  ] 


* « -  •  -  J  ■*.:  <.'.j  ;  j<i.  h  . .  "*  ?»C.O  :  j/.  ’  !^'1.  : SSC.D 


Figure  4.1a.  Location  of  laser  lines  2v[R(12)],  2v(R(l4)],  ?v(R(16)] 
and  ?v[R(lfl)]  in  the  spectral  region  under  investigation. 
Table  4.1  describes  the  experimental  conditions  under 
which  each  row  of  spectra  was  recorded. 
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Figure  4.1b.  Location  of  the  laser  line  2v(R(20)]  in  the  spectral 
region  under  investigation.  Table  4.1  describes  the 
experimental  conditions  under  which  each  row  of  spectra 
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was  recorded 


to  be  five  times  [28]  the  AFGL  line  width  listing,  then  the  2v[R(12)]  is 
1.38  line  widths  away  from  line  center.  The  2v[R(14)]  =  1943.8618  cm*l 
line  is  1.3448  cm“l  away  from  a  very  strong  line  centered  at 
1942.5170  cm-1.  Using  the  same  assumption  as  above,  the  2v[R(14)]  line 
is  estimated  to  be  5.75  line  widths  from  line  center. 

As  for  the  2v[R(18)]  =  1949.2452  crW  line,  there  is  no  doubt  that 
it  is  well  within  the  line  width  of  a  local  line  centered  at 
1949.2500  cm~l.  In  fact,  it  is  0.014  line  widths  away  from  line  center 
using  the  above  line  width  criterion  of  5  times  the  AFGL  line  width 
listing.  This  line  should  give  an  opportunity  to  observe  the  near  line 
center  shape  as  a  function  of  temperature  and  pressure. 

Another  laser  line  in  this  region,  the  2v[R(20)]  =  1951.8622  cm-* 
shown  in  Figure  4.1b,  is  also  of  interest  because  it  could  possibly  be 
paired  with  the  2v[R(18)]  line  as  an  on-line/off-line 
(2v[R(18)]/2v[R(20)])  resonant  pair  for  use  in  differential  lidar  remote 
sensing  of  atmospheric  water  vapor.  At  least  one  frequency  doubled  C0j> 
DIAL  system  is  currently  being  developed  [21]  and  the  characteristics  of 
this  pair  may  be  of  interest  and  thus  should  be  investigated  further. 

In  addition,  these  two  lines  appear  to  be  off  line  resonance  of  the  CO? 
molecule  which  can  often  be  a  problem  in  the  near  infrared  region. 

Observations  of  these  lines  were  made  for  three  temperatures, 

295’K,  310°K,  and  323°K.  The  323’’K  temperature  is  as  high  as  the  safe 
operating  upper  limit  of  the  absorption  cell  insulation.  Temperatures 
lower  than  room  temperature  are  also  of  interest  but  this  study  chose 


to  focus  on  room  temperature  and  above. 

The  water  vapor  partial  pressures  were  generally  5,  10  and  15  torr. 
Several  data  points  at  295°K  were  taken  with  water  vapor  pressures 
slightly  different  than  the  three  listed  above  due  to  an  error  in 
selecting  the  proper  volume  of  H2O  samples  to  evaporate  into  the  cell. 
The  error  was  detected  after  the  absorption  measurements  were  taken  and 
the  actual  pressures  had  to  be  recalculated  and  are  believed  to  be 
correct  for  those  particular  measurements  at  295°K. 

C.  DATA  PRESENTATION  AND  ANALYSIS 

1 .  Data  Presentation 

Tables  4.2  through  4.4  present  the  experimentally  measured 
absorption  coefficients  as  a  function  of  water  vapor  pressure  for  the 
investigated  laser  lines  and  temperatures .  Along  with  the  data  are  the 
theoretical  absorption  coefficient  values  determined  by  use  of  the  TN 
model.  Both  the  experimental  and  theoretical  values  are  plotted  for 
each  laser  line  and  temperature  in  Figures  4.2  through  4.10. 

The  error  bars  in  Figures  4.2  through  4.10  result  from  a  shift  in 
the  measured  background  (or  100%  transmittance  level)  during  the 
experiment.  One  of  the  end  points  of  each  error  bar  represents  the 
absorption  coefficient  value  based  on  the  intitial  background 
measurement  and  the  other  is  based  on  the  final  background  measurement, 
after  the  sample  has  been  evacuated  from  the  cell.  As  discussed  in 


TABLE  4.2 


NITROGEN  BROADENED  H20 

ABSORPTION 

DATA  GATHERED  FOR  THE 

2v(R(12) )  =  1941 

.096  cm“l 

:02  LASER 

LINE 

h2o* 

Theoreti cal 

Experimental 

Temperature 

Pa 

k (km-^ ) 

k  (k 

m-1 ) 

(°k) 

5.4 

5.0546 

5.5680  - 

5.8765 

295 

8.1 

7.7285 

7.7529  - 

8.3540 

295 

10.0 

9.6686 

10.1480  - 

10.8090 

295 

16.1 

16.2226 

16.1225  - 

16.7310 

295 

5.0 

4.7765 

5.3859  - 

5.4317 

310 

10.0 

9.8850 

10.8751  - 

11.4873 

310 

15.0 

15.3243 

17.0015  - 

17.2309 

310 

5.0 

4.8561 

_ 

_ 

323 

10.0 

10.0411 

7.8369  - 

9.0587 

323 

15.0 

15.5540 

“ 

- 

323 

„ . 

*A1 1  samples  nitrogen  broadened  to  740  Torr  total  pressure. 
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TABLE  4.3 


NITROGEN  BROADENED  H20  ABSORPTION  DATA  GATHERED  FOR  THE 
2v(R(14))  =  1943.862  cm'1  C02  LASER  LINE 


h2o* 

Theoretical  Experimental 

Temperature 

Pa 

k (km-1) 

k(k 

m-i) 

(°k) 

5.4 

5.1312 

6.1589  - 

6.2922 

295 

8.1 

7.8440 

8.6566  - 

9.5165 

295 

10.0 

9.8117 

11.1056  - 

11.5271 

295 

16.1 

16.4561 

17.0808  - 

17.9375 

295 

5.0 

4.6333 

5.4146  - 

5.7333 

310 

10.0 

9.5895 

11.3992  - 

12.2251 

310 

15.0 

15.3243 

16.7991  - 

17.1013 

310 

5.0 

4.5324 

4.2527  - 

4.6338 

323 

10.0 

9.3760 

7.1591  - 

7.9797 

323 

15.0 

14.5298 

- 

- 

323 

*A11  samples 

nitrogen 

broadened  to 

740  Torr 

total  pressure. 
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TABLE  4.4 


NITROGEN  BROADENED  H20  ABSORPTION  DATA  GATHERED  FOR  THE 
2v(R(18))  =  1949.245  cm-1  C02  LASER  LINE 


h20* 

Theoretical 

Experimental 

Temperature 

Pa 

k (km-^ ) 

k(km_l ) 

(°k) 

5.4 

6.0398 

6.5052  -  7.0831 

295 

8.1 

9.8044 

10.1327  -  10.8879 

295 

10.0 

12.1183 

12.0633  -  12.6428 

295 

16.1 

19.5762 

18.4152  -  19.2715 

295 

5.0 

7.5425 

8.1124  -  8.6830 

310 

15.1128 

15.8347  -  16.8073 

15.0 

22.7013 

23.3303  -  23.6706 

310 

8.9802 

_  „ 

323 

18.0330 

16.5227  -  17.1635 

323 

15.0 

27.0681 

- 

323 

*A ! 1  samples  nitrogen  broadened  to  740  Torr  total  pressure 
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ABSORPTION  COEFFICIENT  (km 


30 


1941.0958cm" 


- 

EXPERIMENTAL 

a  THEORETICAL 

- 

T  =  310  K 

I 

— 

X 

- 

A 

1 
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1 

A 

_ I _ 
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6  8  10  12  14  16 

H20  PARTIAL  PRESSURE  (  Forr) 


■1.3  Comparison  of  expor imonta 1  and  theoretical  'IpO  partial 
pressure  dependence  for  2f[R(I2)]  =  1941. 095,9  c.m~^  at. 
310  K'’. 
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1941.0958  cm 


Figure  4.4  Comparison  of  experimental  and  theoretical  H^O  partial 
pressure  dependence  for  2v[R(12)]  =  1941.0958  cm~l  at 


1943.8618cm 


4.S  Comparison  of  experimental  and  theorptical  H^O  partial 
pressure  dependence  for  2v[R(14)]  =  1943.8618  cm-1  at 


ABSORPTION  COEFFICIENT  (Km' 


3 


1949.2452  cm 


20 


10 


1  EXPERIMENTAL 
A  THEORETICAL 
T  =  295  K 


I 

A 


1 


4  6  8  10  12  14 

H20  PARTIAL  PRESSURE  (Torr) 


16 


are  4.8  Comparison  of  experimental  and  theoretical  830  partial 
pressure  dependence  for  2v[R(18)]  =  1949.2452  cm~l  at 
295  K°. 
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ABSORPTION  COEFFICIENT  (km" 


H20  PARTIAL  PRESSURE  (Torr) 


Figure  4.10  Comparison  of  experimental  and  theoretical  H2O  partial 
pressure  dependence  for  2v[R(18)]  =  1949.2452  cm-*  at 
323  K° . 
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Chapter  III,  establishing  a  stable  background  was  the  major  and 
virtually  unresolved  problem  encountered  throughout  the  study.  The 
problem  finally  got  so  dominant  that  only  a  small  amount  of  data  was 
able  to  be  obtained  at  323°K.  The  data  that  was  obtained  at  the 
2v[R(12)]  and  2v[R(18)]  lines  is  questionable  and  cannot  be  taken  with 
much  confidence.  Data  measured  for  295°K,  310°K  and  the  2v[R(14)]  line 
at  323°K  is  believed  to  be  more  accurate  but  the  instabilities  of  the 
system  greatly  de-sensiti zed  all  of  the  absorption  measurements. 

2.  The  Empirical  Expression  for  the  Absorption  Coefficient 

The  general  form  of  the  absorption  coefficient  which  is  used  to 
express  the  data  is 


,)  = 

RT 


(pn2cn 


+  paCS^  ^km 


(42) 


where  C$  is  the  self-broadening  coefficient  for  water  vapor,  C^  is  the 
nitrogen  broadening  coefficient,  p^  i s  the  pressure  of  nitrogen  in 
torr,  p  is  the  water  vapor  partial  pressure  in  torr,  R  is  the  ideal 

a 

gas  constant,  and  T  is  the  temperature  in  Kelvin.  For  this  data,  the 
following  expression  is  also  true: 


pT  =  PN?  +  Pa  =  740  torr  . 


(43) 
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(44) 


Substituting  for  in  Equation  (42)  yields, 

M»)  -Hf  ((PT  -  pa)  cN  -  Pacs)  [km-1] 

k(v)  -  (£TCf|)pa  +  {£S^N)p|  [km-1]  .  («) 

Thus,  a  general  second  order  polynomial  can  represent  the  data, 

k ( v)  =  Apa  +  Bpa2  [km-1]  ,  (46) 


where 


A  =  [km-*  torr"*] 

RT 

8  =  ^S~Cn  [km-*  torr-^] 

RT 

This  expression  is  valid  for  absorption  coefficients  at  frequencies 
further  than  one  line  width  from  a  local  line  center  such  as  appears  to 
be  the  case  for  the  2v[R(12)]  and  the  2v[R(14)]  lines.  However,  the 
self  broadening  term  of  Equation  (42)  becomes  negative  for  frequencies 
within  one  line  width  of  a  local  line  center.  This  leads  tc  a  general 
second  order  polynomial  of  the  form, 

k(v)  =  Apa  -  Rpa2  [km-*]  ,  (47) 

where 

A  =  PJ£n  [km-*  torr-*] 

RT 

8  =  [km-*  torr"^] 

RT 
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The  2v[R(18)]  line  should  follow  the  form  shown  in  Equation  (47)  since 
it  is  within  the  line  width  of  a  local  absorption  line. 

Thus,  the  absorption  coefficient  as  a  function  of  water  vapor 
pressure  for  a  frequency  within  a  line  width  of  line  center  will  have 
linear  or  parabolic  down  characteristics  (negative  second  order 
coefficient);  whereas  the  absorption  function  for  frequencies  outside 
of  a  line  width  will  display  a  linear  or  parabolic  up  characteristic 
with  the  curve  becoming  more  parabolic  as  C$  becomes  much  larger  than 
C|\|.  Equation  (42)  can  be  expressed  yet  another  way, 

k(v)  =£n  pa  (pN  +  8pa)  [km-1]  (48) 

where 


B  is  large  in  the  window  regions  where  absorption  is  due  to  the  wings  of 
far  away  lines.  In  the  10  gm  region,  for  example,  B  ranges  from  200  to 
300.  The  quadratic  term,  therefore,  dominates  in  this  region  and  the 
absorption  coefficient  versus  water  vapor  partial  pressure  is  parabolic. 

In  the  high  absorption  region  of  this  study,  where  the  absorption 
is  due  primarily  to  local  lines,  B  is  small  (5  -  10)  and  thus  the 
absorption  coefficient  versus  water  vapor  partial  pressure  functions  are 
nearly  linear. 
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.  A  Least  Squared  Curve  Fit  to  the  Data  and  Analysis 


Using  the  general  second  order  polynomial  form  of  the  absorption 
coefficient  given  in  Equation  (46),  a  least  squared  curve  fit  was  made 
for  the  experimental  data  displayed  in  Figures  4.2  through  4.10.  The 
resulting  equations  are  presented  in  Table  4.5  along  with  the  least 
squared  fit  equations  determined  from  the  TN  model  theoretical  values 
given  in  Figures  4.2  through  4.10. 

The  absorption  coefficient  functions  of  this  heavy  local  line 
absorption  region  appear  nearly  linear  as  predicted  by  the  TN  model  and 
observed  by  past  studies  of  heavy  absorption  regions.  The  theoretical 
second  order  term  is  two  orders  of  magnitude  smaller  than  the  linear 
term  for  the  2v[R(12)]  =  1941.096  cm_l  line  and  the  2v[R(14)]  = 

1943.862  cm~l  line,  and  four  orders  of  magnitude  smaller  for  the 
2v[R(l8)]  line  which  is  much  closer  to  line  center.  The  data  appears  to 
agree  with  the  theoretical  functions  in  at  least  the  linear  term  for  the 
2v[R(12)]  and  the  2v[R(18)]  lines  at  295°K  and  310°K.  The  agreement  in 
the  linear  term  with  the  model  for  the  2v[R(14)]  data  seems  to  be  less 
close  than  for  the  other  two  lines.  All  of  the  linear  terms  for  the 
data  at  295°K  and  310°K  are  consistently  higher  than  the  TN  model 
predicts.  Since  there  is  such  a  small  amount  of  accurate  data  to  go  on, 
it  is  difficult  to  say  where  the  inaccuracy  lies  -in  the  data  or  in  the 
model . 

As  for  the  quadratic  terms  of  the  experimental  measurements,  they 
seem  to  differ  in  sign  from  their  theoretical  counterparts  for  all 
except  the  2v[R(12)]  line  at  310°K.  The  difference  between  the 
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TABLE  4.5 


EXPERIMENTAL  AND  THEORETICAL  ABSORPTION  COEFFICIENT 
k(km-l)  =  Ap  +  Bp2  WITH  WATER  PRESSURE  p  IN  TORR  FOR 
THE  INVESTIGATED  LASER  LINES 


EXPERIMENTAL  THEORETICAL 

WAVENUMBER  TEMPERATURE  ABSORPTION  COEFFICIENT  ABSORPTION  COEFFICIENT 
(cm-1)  (°K)  EQUATION  EQUATION 


1941.096 

295 

1941.096 

310 

1941.096 

323 

1943.862 

295 

1943.862 

310 

1943.862 

323 

1949.245 

295 

1949.245 

310 

1949.245 

323 

(1.05)p-(l. 95E -3 '  2 
(1.09)p+(3.24E-3)p2 
(8.53E-1 )p+(0.0)p2 

(1.19) p-(6.68E-3)p2 

(1.19) p-(3.63E-3)p2 
(1.09)p-(3.30E-3)p2 

( 1 • 42 )p - ( 1 .62E-2)p2 
(1.72)p-(9.95E-3)p2 

(1 ,68)p+(0.0)p2 


(9.00E-1 )p+(6.68E-3)p2 
(9.92E-1 )p+(5.63E-3)p2 
(9.33E-1 )p+(6.91E-3)p2 

(9.14E-l)p+(6.72E-3)p2 
(8.94E-1 )p+(6.45E-3)p2 
(8.75E-l)p+(6.2IE-3)p2 

(1.20)p+(7.06E-4)p2 
(1.51 )p+(4.61E-4  )p2 
(1.80)p+(5.97E-4)p2 
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experimental  and  theoretical  quadratic  terms  may  he  within  the 
experimental  error  of  the  measurements.  The  instability  of  the 
background  measurement  has  made  the  absorption  measurement. s  too 
imprecise  to  determine  this  small  quadratic  term  with  confidence. 

The  data  measured  at  323°K  illustrates  the  di f f ( cu 1 *• i es  encountered 
toward  the  end  of  the  experiment.  The  background  measurements  oecame 
less  stable.  As  a  result,  not  enough  data  was  qathered  to  even  estimate 
a  quadratic  term  for  the  2  v[R  (12)1  and  2vt.R(18)1  lines.  As  Figure  4.8 
indicates,  the  most  data  gathered  for  this  temperature  was  for  the 
2v[R(14)]  line  and  indeed  Table  4.5  shows  this  data  pretty  much  falling 
into  the  same  range  of  agreement  as  data  at  the  lower  temperatures .  In 
fact,  the  data  for  this  2v[R(I4)]  line  at  323°<  with  respect  to  its 
linear  term,  seems  to  agree  with  the  theoretical  linear  term  even  better 
than  2v[R(14)]  data  measured  at  the  lower  temperatures. 

Although  the  data  has  a  wide  margin  of  uncertainty,  the  overall 
indications  seem  to  be  that  it  is  substantiating  the  TN  model  line 
shape.  Flowever,  it  is  clear  that  much  more  sensitive  data  must  be  taken 
before  any  definite  improvements  can  be  made  in  the  model  resulting  from 
observations  of  the  local  line  shapes  as  a  function  of  temperature  for 
this  region.  Until  such  time  as  more  accurate  data  is  obtained,  it  is 
pointless  to  analyze  this  data  any  further. 


78 


CHAPTER  V 


CONCLUSIONS 

Water  vapor  absorption  in  the  window  regions  of  the  near  infrared 
is  of  considerable  interest  to  many  infrared  energy  transmission 
systems.  A  controversy  exists  as  to  what  the  cause  or  causes  are  for 
the  water  vapor  continuum  absorption  observed  in  these  window  regions. 
One  school  of  thought  says  that  it  is  a  result  of  complex  clusters  of 
water  molecules  (i.e,  dimers)  that  absorb  radiation  as  a  slow  varying 
function  of  frequency  throughout  the  entire  near  infrared  region. 
Important  to  the  support  of  this  dimer  theory  is  the  fact  that  it 
explains  the  negative  temperature  dependence  of  the  continuum. 

Another  school  of  thought  says  that  the  continuum  is  from  the  far 
wing  overlap  into  these  window  regions  of  very  strong  absorption  lines 
centered  deep  in  the  vibrational  bands  of  H2O.  A  line  shape  model  which 
is  Lorentzian  out  to  5  cm"*  from  line  center  and  modified  in  the  far 
wings  was  developed  at  Ohio  State  University  in  1976-79  by  Michael  E. 
Thomas  and  Dr.  Robert  J.  Nordstrom.  The  form  of  the  far  wing  shape  was 
developed  theoretically  and  then  certain  line  shape  parameters  were 
derived  by  trial  and  error  fit  to  data  taken  in  regions  of  high 
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absorption  dominated  by  local  absorption  lines.  Thus,  the  model  was 
developed  in  regions  where  the  continuum  is  of  little  significance  and 
therefore,  the  development  was  independent  of  any  continuum  effects.  It 
was  then  applied  to  the  window  regions  and  was  the  first,  line  shape 
model  to  successfully  model  the  pressure  and  temperature  characteristics 
of  the  continuum.  Although  the  model  fell  short  of  modeling  the 
magnitude  of  the  negative  temperature  dependence,  it  did  show  a  negative 
dependence  in  temperature  for  the  continuum  absorption. 

The  purpose  of  this  study  was  to  gather  temperature  dependence  data 
in  a  region  dominated  by  local  line  absorption  so  that  the  line  shape 
parameters  used  in  the  Thomas-Nordstrom  model  can  be  improved  on  and 
more  insight  can  be  made  as  to  the  temperature  dependence  of  these  line 
shapes.  The  study  uses  a  frequency  doubled  C0;>  laser  frequency  probe  to 
examine  the  5.15  urn  infrared  region  at  three  frequency  doubled  lines: 

1)  2v[R ( 12 ) ]  =  1941.0952  cm"1 

2)  2 v[R ( 14 ) ]  =  1943.8624  cur1 

3)  2v[R ( 18) ]  =  1949.2495  cm"1  . 

These  lines  were  examined  at  temperatures  of  295°K,  310°K,  and  323°K 
with  the  sample  pressures  set  to  5,  10,  and  15  torr.  All  of  the  samples 
were  nitrogen  broadened  to  740  torr. 

Unfortunately,  no  solid  conclusions  or  new  insights  can  be  made 
from  the  data  gathered  during  this  study  due  to  the  experimental 
instabilities  of  the  optical  and  detection  systems  and  therefore,  the 
inaccuracy  of  the  measurements.  It  would  appear  that  the  data  supports 
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the  gross  features  of  the  models'  ability  to  predict  the  absorption 
coefficient  versus  water  vapor  partial  pressure  for  the  temperatures 
investigated.  It  is  clear,  however,  that  much  more  accurate  data  must 
be  measured  for  these  lines  before  any  real  improvements  in  the  model 
can  be  made.  The  instability  of  the  system  and  thus  the  inability  to 
maintain  a  constant  background  throughout  the  measurements  was  the  major 
problem  encountered  in  this  first  attempt  to  obtain  this  data. 

These  instabilities  were  identified  in  two  areas.  First  it  was 
discovered  that  the  signal  detector  was  giving  fluctuating  readings 
seemingly  dependent  on  how  much  nitrogen  remained  in  its  dewar.  These 
fluctuations  were  independent  of  the  COj  laser  and  the  second  harmonic 
generation  because  the  laser  source  during  these  observations  was  a 
Helium  Neon  laser.  When  the  signal  detector  was  replaced  by  a  HgCdTe 
detector,  which  was  the  only  other  detector  available,  a  stable 
background  was  measured.  It  is  highly  recommended  that  the  signal 
detector  be  replaced  with  a  new  InSb  detector. 

The  other  instability  identified  was  in  the  CO^  laser.  Using  the 
HgCdTe  detector  as  the  signal  detector,  only  the  CO?  laser  was 
reintroduced  into  the  system  and  the  background  measurements  were  again 
unstable.  It  was  never  determined  how  this  laser  could  produce 
fluctuations  in  the  background,  but  it  is  recommended  that  the  Sylvania 
CO?  laser  previously  described  [26]  be  used  as  the  radiation  source. 

The  Sylvania  laser  has  a  cw  power  output  of  2  to  3  watts  which  has  been 


found  to  be  sufficient  to  produce  a  measurable  amount  of  second 
harmonics.  At  the  beginning  of  the  project,  it  was  believed  that  10 
watts  cw  would  be  required  in  order  to  produce  a  reasonable  amount  of 
second  harmonics  and  thus  the  reason  for  constructing  the  new  laser.  As 
it  turned  out,  in  order  to  be  assured  of  a  TEMqq  mode  for  the  output  of 
the  new  laser,  an  iris  was  required  in  the  cavity  which  decreased  the 
laser's  power  to  the  same  range  as  the  Sylvania  laser.  The  Sylvania  is 
the  better  of  the  two  lasers  for  this  output  power  range  and  should 
therefore  be  brought  back  into  working  order  and  put  into  the  optical 
system. 

The  exact  role  in  the  instability  problem  of  the  CdfieAs?  crystal 
and  nitrogen  dewar  was  never  determined.  It  was  found  that  the  dewar 
had  to  be  suspended  in  order  to  keep  the  outer  skin  thermal  expansion 
and  contractions  of  the  dewar  from  disorienting  the  crystal.  This 
problem  was  encountered  when  the  dewar  was  resting  on  a  pedestal  and 
supported  by  its  tail  piece. 

The  crystal's  face  has  been  pitted  by  exposure  to  radiation  of  too 
much  intensity  and  there  only  seems  to  be  one  spot,  on  the  crystal  face 
where  a  fair  amount  of  second  harmonic  energy  is  observed  when  exposed 
to  the  CO 2  radiation.  Therefore,  getting  a  new  crystal  would  be  the 
best  direction  to  take  in  constructing  a  better  system.  Now  that  the 
crystal  focusing  optics  have  been  developed  there  would  be  no  need  to 
expose  a  new  crystal  to  high  input,  power  intensities  and  so  crystals 
used  in  the  future  should  be  safe  from  face  surface  damage. 
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Once  these  improvements  in  the  system  are  made,  data  should  be 
retaken  for  the  three  laser  lines  investigated  in  this  study  at  similar, 
if  not  the  same,  temperatures.  If  this  data  appears  to  be  accurate 
enough,  the  C$  sel f-broadeni ng  coefficients  should  be  determined  from 
the  data  and  their  characteristics  as  a  function  of  temperature 
observed.  In  addition,  the  2v[R(20)]  line  should  be  investigated  as  a 
possible  partner  in  an  on-line/off-line  resonance  differential 
absorption  HOAR  pair  with  the  2v[R(lS)l  line.  The  2v[R(18)]  laser  line 
is  well  within  the  halfwidth  of  a  fairly  strong  local  absorption  line 
and  the  2v[R(20)]  line  sits  within  a  narrow  window  of  this  5.15  ym 
region. 
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